NONLINEAR WAVE EQUATIONS AND SINGULAR SOLUTIONS 



HIDESHI YAMANE 

Abstract. We construct solutions to nonlinear wave equations that are sin- 
gular along a prescribed noncharacteristic hypersurface which is the graph of 
a function satisfying not the Eikonal but another partial differential equation 
of the first order. The method of Fuchsian reduction is employed. 



1. Introduction 

In our previous article |10j . we constructed a solution to a nonlinear Cauchy 
problem and gave a lower bound of its lifespan. 

In the present article, we consider certain nonlinear wave equations with poly- 
nomial nonlinearity. More precisely, the right hand side of an equation in our class 
is a polynomial of degree 2 in the derivatives of the unknown function. Notice that 
some of this kind of equations fall within the class treated in |1U) . 

We construct solutions that are singular along a prescribed noncharacteristic 
hypersurface satisfying a certain condition written in terms of a partial differential 
equation of the first order. Blowup occurs due to the presence of a logarithmic 
term. Its coefficient (and hence the speed of blowup) is determined by the hyper- 
surface and the nonlinear term of the equation. On the other hand, it is possible to 
construct hypersurfaces in accordance with a given pair of the coefficient and the 
nonlinear term. 

In later sections, we briefly discuss a nonlinear elliptic equation and multi-valued 
solutions to partial differential equations in the complex domain. We also consider 
equations with nonlinearity of higher order. 

The proofs are based on Fuchsian reduction. This method, a powerful tool in 
many branches of analysis, was first employed by Kichenassamy-Littman in their 
construction of singular solutions to the equation Ou = exp u in 0] . Similar meth- 
ods were used in 0, and [Hj- We reduce our equations to Fuchsian ones and 
use results by Gerard- Tahara in pQ and (Chapter 8) . An alternative proof is the 
reduction to first-order Fuchsian systems which are solved in 0] and 

2. Main result 

Let f(t, x; t, £) be a polynomial of degree 2 in (r, ()elx R™ whose coefficients 
are analytic functions in (t,x) in an open set U of R x R™. We decompose / 
into homogeneous parts: / = /a + /i + fo, where // is homogeneous of degree /. 
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For I = 1,2, we sometimes denote fi(t,x;T,£) by fi(r,£) for brevity. Notice that 
/o = /o(t) x) is just a function in (i, x) and is independent of (r, £). 

The primary purpose of the present paper is the study of a nonlinear wave 
equation 

2 

(2.1) u tt -Au= f(t,x;d t u,Vu) = y^f l (t, x;d t u,Vu), 

1=0 

where V is the gradient with respect to x. 

We construct a solution u(t, x) which blows up on a prescribed noncharacteristic 
hypersurface. More precisely, our main result is the following: 

Theorem 2.1. Let E be a hypersurface defined by E = {(t,x) £ U;t = ip(x)}, 
where ip(x) is an analytic function. Assume 

(2.2) 0^1- \V^(x)\ 2 = af2(ip(x),x; -1, 

in U for a nonzero constant a. Then for any analytic function vq onT,, there exists 
an analytic function v(t, x) in an open neighborhood Q of E such that: 

• u(t, x) — —a log(t — ip(x)) + v(t, x) is a solution to $2.1\) in fl n {t > ip(x)}, 

• the restriction of v on £ coincides with vq . 

Such a function v(t,x) is unique in a neighborhood o/E. 

Remark 2.2. [A] If / 2 (t,£) = a^r 2 - |£| 2 ), then |22J) reduces to l-|V?A(a;)| 2 ^ 0. 
It is satisfied by almost any function ip. 

[B] The function if) is a solution not to the Eikonal equation but to another equa- 
tion of the first order. It would be an interesting task to formulate a general theory 
encompassing the "pseudo-Eikonal" equations (|2.2() . I|8.2|l and the curvature condi- 
tion in 0j. A promising approach is to deal with them in the framework of but 
it is beyond the scope of the present paper. 

[C] The speed of blowup is controlled by the coefficient a. It is calculated in terms 
of /2 and ip. On the other hand, we can find many ^'s satisfying 1)2. 2|) for a given 
pair (/2, a). We have only to use the Cauchy-Kovalevskaya theorem or the classical 
theory of partial differential equations of the first order. 

Recall that the speed of blowup was always the same irrespective of ij) in the 
study of Du = e u in g] and [H] • 

3. Proof of Theorem 12. II 

We introduce a new system of coordinates (T, X) by T = t — ip(x), Xi = Xi(i = 
1,2,..., n). The derivatives in the new variables are marked with" : 

d t = d/dx h d t = d/dt, d t = d/dXt, d x = 0!, . . . , d n ), d T = d/dr. 

In the same way, A and V are the Laplacian and the gradient in x respectively, 
while A and V are those in X. Subscripts on functions refer to the derivatives in 
x. For example, ipi = diip,ipu = dfip- 

The derivatives transform according to the following rules: 

d t =d T , di = -ipid T + c\. 

We want to express the d'Alembertian □ = df — A in terms of di and dr- This 
calculation has already been explained in 0] . We nevertheless repeat the argument 
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here for the sake of the reader. It might seem a little confusing since it employs 
two systems of coordinates simultaneously. To begin with, we get 

Of = di(-ipA + di) = -di^A) + dA 

= -ijjididr - 4>udr + (-ipidr + di)di 

= -ipi(-4>idr + di)d T - ipud T - ^d T di + df 

= -0 2 <9| - 2ipidid T - ipud T + df- 

Therefore □ has the following form: 

n 

(3.1) □ = + 2 + (A^)9 T - A, 

i=l 

where * = f(X) = $(x) = 1 - \Vip(x)\ 2 ^ 0. 

Set u(t,x) = — alogX 1 + v, where v is an unknown analytic function near £ = 
{T = 0}. Then by the left hand side of JTTJ) is 

(3-2) □„=___£+□„. 
Next, let us calculate the right hand side. We have 

dtu = —a/T + drv, diU = aipi/T — ipidrv + diV. 



It leads to 



(d t u) 2 = ^~^d T v + (d T v) 2 , 



(d t u)(d t u) = + "^j^Btv - ?-diV + (d T v)(-ipAv + d.v), 

{diu)(dju) = a 'pf 1 + ^ {-2aipiipjd T v + a(ipidj + ipjdi)v} 

+ (-ipAv + d l v)(-tp j d T v + djv). 

On the other hand, there exists an analytic function f 2 {T, X) in a neighborhood of 
£ such that 

f 2 (t,x;-l,ViP(x)) = f 2 (iP(x),x;-l,ViP(x)) +Tf 2 (T,X). 

We set f 2 (T,) = f 2 (ip(x),x; — l,Vip(x)) for brevity. This quantity appears in (|2.2|) 
and should be distinguished from f 2 (— 1, Vip(x)) — f 2 (t, x; —1, \7ip(x)), a particular 
case of fi{r, £) = f t (t, x; r, £)■ We get 

M-l^(x)) _ / 2 (S) h(T,X) 

rp 2 T 2 T 

Consequently, the right hand side of the equation 1|2.1|) is 

(**\ tU X V7 \ « 2 / 2 (S) , a 2 .f 2 (T,X) 2a ~ 
(3.3) f{t, x; d t u, Vu) = — — — + - j^f2{-l, v V0<9tv 

+ ^L(T, X, d x )v + h (d T v, {-i/h&tv + 

+ (d T V,{-lPAv + diV}^+f (t,x), 

where L(T, X, dx) is a vector field with [L, T] = 0. 
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Now let us compare ()3.2() and ((3.3(1 . Our assumption ((2.2(1 implies that the 
coefficients of T~ 2 coincide: a\t = a 2 /2(S). Therefore, in view of JH3J, the equation 
(|2.1|) is equivalent to the following one: 

(3.4) - aAiP + T | + 2 ^ ^a^T^ + (&i>)d T v - Auj 

= a 2 / 2 (T, X) - 2*cVf + aL(T, X, d x )v + Tf 2 (d T v, {-^d T v + 8,v}?j 
+ a/i(-l,V^)+r/i (ar^{-^9T« + 4«}i) +Tf (t,x). 

We divide this equation by "J ^ and get an example of (|4.1|) in 21 below. By 
Proposition ^. II there exists a unique analytic solution v to l|3.4|l with u| s = vq in 
a neighborhood of E. The function u = — alogT + v solves (I2.1|) in {T > 0}. We 
have finally proved Theorem 12. II 



4. Nonlinear Fuchsian partial differential equations 

In this section, we prove a result on nonlinear Fuchsian partial differential equa- 
tions which was used in the proof of Theorem 12. II 

We consider a nonlinear equation in a neighborhood of {0} x V C R* x K™, 
where V is an open set of R™. We employ the usual notation for derivatives: 
dt = d/dxi, d x = (dt, ■ . . , d n ),Vv = (dtv, d n v) and d t = d/dt. 

Now let us consider a Cauchy problem for the following nonlinear Fuchsian partial 
differential equation of weight 1 : 

(4.1) td 2 v + 2d t v = a(t, x) + tP(t, x, d x )d t v + Q(t, x, d x )v 

+ tj3(t, x)(d t v) 2 + R{t, x, d x )v ■ td t v + S(t, x; Vv). 

Here a and (3 are analytic functions. The linear differential operators P, Q and R 
with analytic coefficients are of order 1,2 and 1 respectively and [P,t] = [Q,t] = 
[R, t] = 0. We assume that S(t, x; £) is a quadratic form in £ whose coefficients are 
analytic functions in (t, x). 

Proposition 4.1. Let vq(x) be an analytic function in V. Then (|4.1() has a unique 
analytic solution v(t, x) with v(0, x) — vq(x) near {0} xl^Ctx R™. 

Proof. First we reduce the problem to the case vq(x) = 0. We introduce a new 
unknown function w by v(t,x) = vq(x) + w(t,x). Then we have td 2 v + 2dtV = 
td 2 w + 2d t w. 

We denote the right hand side of 1(4.1(1 by N[v]. An elementary computation 
shows 

N[v] = N[v + w] 

= N[w] + Qv + Rv ■ td t w + S(t, x; Vv ) + T(t, x; Vw), 

where T(t, x;£) — 2<=i c *(^ x )£i f° r some analytic function Ci(t,x). This implies 
that w(t,x) satisfies an equation of the type 14.1(1 . possibly with different coeffi- 
cients. Therefore we have only to prove that 14.1(1 has a unique analytic solution 
v(t,x) with v(0,x) = 0. Multiplication by t yields a nonlinear Fuchsian partial 
differential equation solved in pQ and [2]. Its characteristic polynomial is p 2 + p 
and there exists a unique analytic solution v(t, x) with v(0, x) = 0. □ 
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Remark 4.2. [A] There is a stronger uniqueness result due to Gerard- Tahara. We 
only give a sketch here and refer the reader to ^ or [2] for a precise statement. 

Assume that vo(x) = for simplicity and let v(t, x) be a multi- valued holomor- 
phic solution to (|4.1() in the universal covering space of {t ^ 0} C Cf x C™. If 
there exists a positive constant p such that v(t, x) = 0(|t| p ) as t tends to zero, then 
v(t,x) is nothing but the solution v(t, x) in Proposition l4.il 

This observation enables us to improve the uniqueness part of Theorem 12.11 It 
rules out some functions involving terms of the form T a log' 3 T, a > 0, T — t — ip(x). 
[B] The present formulation, a very restrictive one, has been chosen because it is 
all we need in the proof of Theorem 12.11 There is much room for generalization. 
We can easily consider other characteristic polynomials than p 2 + p. Difficulties 
arise if we consider nonlinear terms which are not entire in its arguments; it is not 
always possible to consider large initial values. 

5. Time- reversal 

In this section we discuss time-reversal and construct solutions on the negative 
side of £ under certain conditions. The equation ((2.1(1 is not invariant under time- 
reversal. Indeed, if we set s — —t, it becomes 

(5.1) u ss - Au = f(—s, x. ~d s u, Vit). 
In addition to (|2.2Jl . we assume 

(5.2) f 2 (t,x;l,W) =/ 2 (*,z;-l,V# 

Here notice that /2(t, £) = fi{— r, — £) is trivial. The condition 15.211 is satisfied if 
/2C?") is free of terms of the form cj(t, x)r£j. 

We set ip(x) — —i/j(x) and introduce /2(r, £) = /2(— t, £), which is, to 1)5.1(1 . what 
/2 is to (12.1(1 . Combination of 1(2.2(1 and 1(5.2(1 yields 

(5.3) 0^1- |V^(x)| 2 = o/ a (-l, V^(x)). 

By using Theorem 12. II we can find solutions to ((5.1(1 in {(s, x); s > ip(x)}. It gives 
a solution to ((2.1(1 in {(t, x); i < V'( :r )}i the negative side of E. 

To sum up, under 1(2.2(1 and 1(5.2(1 . there exist solutions defined on both sides of 
the noncharacteristic hypersurface £ and their singularities propagate along it. 

6. Elliptic case 
We can apply a similar argument to 

(6.1) Au = a" 1 |Vu| 2 

in R™, where a is a nonzero constant. In this case, we can treat any analytic non- 
singular hypersurface A at least locally. Indeed, after rotation, it can be written in 
the form A = {x± = 4>{x2, ■ ■ • , x n )}. Notice that the equation is rotation-invariant. 
The function 1 + J^^i^^/ ® x i) 2 > replaces "J. Hence we have the following: 

Theorem 6.1. Let p be a point of an arbitrary analytic nonsingular hypersurface 
A. Then there exists a neighborhood U of p such that in a connected component of 
U \ A there exists an analytic solution u{x) to \6. 1)) which blows up along AnU . 
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7. Multi-valued solutions in the complex domain 

Our calculus is valid in the complex domain and the solution u(t, x) in Theo- 
rem ^, ll bccomes a multi- valued holomorphic function around the noncharacteristic 
complex hypersurface S c = {(t,x) G C n+1 ;t = ip(x)}. It is the unique solution to 
the following singular Cauchy problem near a point (to,p) G S c ,to — ip(p)'- 

Du — f(t, x, dtu, Vm), 
< u(t ,x) = -alog(t - ip(x)) +v(t ,x), 
d t u(t , x) = -a/ (t - ip(x)) + d t v(t , x). 

For such Cauchy data, singularities propagate along the noncharacteristic hyper- 
surface £ and only along it. Notice that some results in 0] and [H] are open to 
similar interpretations. 

Several authors have constructed multi-valued solutions to nonlinear equations 
around characteristic hypersurfaces. See pQ, [2], [5] and other works by Kichenas- 
samy and by Tahara about Fuchsian problems. They give solutions which are 
singular along the initial surface {t — 0}. Moreover, see [§J and the references cited 
therein for singularities obtained by solving Eikonal equations. 



8. Equations with nonlinearity of higher order 

The leading term of a singular solution can be guessed by solving an ordinary 
differential equation. For example, one finds a solution u(t) = log 2 — 21ogi to 
d 2 u/dt 2 = e* and it was the basis of 0] and 0. On the other hand, the solution 
u{t) = — logt to d 2 u/dt 2 = (du/dt) 2 is the simplest case of our Theorem 12. II 

We are naturally tempted to solve d 2 u/dt 2 = —(du/dt) m+1 for m > 2 and get 
u(t) = C m t( m - 1 '>/ m , where C m = m^ m - r >/ m /{m - 1). This solution vanishes as 
t > tends to 0, but its derivative blows up. In this section, we give a multi- 
dimensional version of this fact. 

Let f(t, x; t, £) be a polynomial of degree m + 1 (m > 2) in (r, £) G R x R™ 
whose coefficients are analytic functions in (t, x) in an open set U of R x R™. We 
decompose / into homogeneous parts: / = X)/=o fh where /; is homogeneous of 
degree I. We sometimes denote fi(t,x;T,£) by fi(r,£) for short. 

We consider 

m+i 

(8.1) u tt - Au = f(t,x;d t u,Vu) = ^ fi(t,x;d t u,Wu) 

1=0 

and construct a solution u(t, x) which vanishes on a prescribed noncharacteristic 
hypersurface. 

Theorem 8.1. Let E be a hypersurface defined by X = {(t, x) G U;t = ip(x)}. 
Assume 

(—m -4- 1 \ m n m 

(8.2) 0^1- |V^(x)| 2 = [ m ^ m \ / m+ i(^(a;),a;;-l,VV'(a:)), 

(8.3) f m (i/>{x),x;-l,V1>(x))=0 
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for a nonzero constant a. Then there exists an open neighborhood of E and an 
analytic solution u(t, x) to \H.1\) in fl n {t > ip(x)} such that: 

u(t, x) ~ a(t — ?p(x)^ m m (vanishing), 

dtu(t,x) ~ a (i — ip(x)) (blow-up) 

to 

as t — ip(x) — > +0. 27ms solution u{t,x) is analytic in ((t — if}) 1 /" 1 , x) and no 
logarithm is involved. 

Remark 8.2. The conditions (|8.2(l and (|8.3(l are equivalent to the existence of such 
analytic functions f m +i(T, X) and f m (T,X) that 
m m_1 

<E2) iWi(M;l,-VV^)) = 7 ——;{l~\^{x)\ 2 )+Tf m+1 (T,X), 

' [—m + l) m a m 

631) / m (*,*;l,-V^(z)) =Tf m (T,X). 

Proof. We introduce a system of coordinates (T, X) as before and seek a solution 
m of the form 

(8.4) u = aT i - m - 1 ^ m + Tv, 

where v is an unknown function analytic in (T x / m , X). It may have fractional-power 
singularities as a function in T. By using the formulas 

d T T = Td T + 1, d\T = Td\ + 2d T 

and H3.1|) . we can express the left hand side of l|8.1(l as follows: 

(8.5) Du=L 1 +a^^-(Ai;)T- 1 / m 

m 

n 

+ (Td| + 2d T ) +2j2MTd T di + di) + (A^)(Td T + 1) - T&}v, 
i=i 

where L 1 = a "^™ ~ ^ ^T-( m+1 >/ m 
to 2 

Next we shall compute the right hand side of (|8.1I) . We get 

f ftu = a 2^iT- 1 /"» + (T$ r i; + tO, 
| diU = -Tpid t u + TdiV. 

Notice that the expression of diU contains dtu. The latter quantity plays an impor- 
tant role in the following calculation. We have 

(8.6) fddtu, Vu) = fi(d t u, {-^Pidtu + Tdiv}i) 

= /j(l,-VV0 (d t u) l + R h 

where Ri is a homogeneous polynomial of degree I in dtu and TdiV and is of degree 
< max{0, 1 — 1}(< to) with respect to d t u. It is of degree > 1 with respect to TdiV 
if I > 1. The "spare" factor T in Tc^v helps us since we are trying to get a Fuchsian 
equation; it offsets the negative power of T in <9 t w. 
We have 

(8.7) (d t u) l = (a— -) T- l / m + l(a— -) ' T-^/ m (Td T v + v) + B u 

y to J y to J 
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where Bi is the remainder of the binomial expansion, that is, 



1-2 /,\ / , x j 



Since / = £Xo /j> we have b y ESI and 
(8.8) f{d t u,Vu) 



m+l 



{L 2 + L 3 + / m+ i(l, + Rl /'C 1 ' -VV)(ft«)', 



where 



1=0 



(_ , \ m-t-i 
TO / 

/ TO - 1 \ m 

L 3 = /m+i(l, -VV»(x)) ■ (m+ 1) a T~ 1 (Tdrv + v). 



We set 



L' 2 = f m+1 (*P(x),x: 1,-V^(x)) U 1 ^) 
then there exists an analytic function L 2 (T, X) such that 

L 2 = L' 2 +T- 1/m L^{T,X). 
The assumption Ij8.2|l implies that L\ = L' 2 and that 



L 3 = _^ . ll+lT-ifTBrv + v). 

TO 

Note that (|8.2|l is written in terms of (—1, Vtp) instead of (1, —Viji). By subtracting 
L\ = L' 2 from (|8.1|l . we come very close to a Fuchsian equation in view of l|8.5|l and 
(Oil . 

We set s = T 1 /™, <9 S = 9/<9s in T > 0. Since toT<9 t = s<9 s , we have 

Tdrv + v = — (sd s v + tow), 
to 

<9 t it = — s _1 {a(TO — 1) + s(sd s v + mi))}. 

TO 

Now we multiply the equation Ou — Li = f — L' 2 by m 2 s m /\E' and obtain a nonlinear 
Fuchsian equation which belongs to the class treated in 1 and Namely it has 
the form: 



f s, X; {(sd s + m)div}i, sd s v + tow, Avj 
-s m ~ 1 L 2 {s, X) — (to + l)(sd s v + mv) + sP 2 (s, X; sd s v + mv)(sd s v + mv) 
+ s m Pz (s, X; sd s v + mv, {div}i 



J2Qi(^ X )s max{m ' l ' 1} {a(m-l) + s{sd s v + mv)} 1 . 



1=0 
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Here L2(s,X) is an analytic function and Pi,Pi and P3 are polynomials in the 
arguments coming after the semicolons. Their coefficients and Qi are analytic 
functions in (s,X) near s = 0. The exponent of s immediately after Qi is not 
m — I but max{m — 1,1} because of (|8.3|) . The characteristic polynomial is (p + 
m)(p + m + 1) and there exists a solution v, which is an analytic function in 

( S ,x) = ((t-^) 1 /"\x). a 

References 

[1] Gerard R. and Tahara H., Solutions holomorphes et singulieres d'equations aux derivees 
partielles singulieres non lineaires, Publ. Res. Inst. Math. Sci., 29(1993), 121-151. 

[2] Gerard R. and Tahara H., "Singular nonlinear partial differential equations", Vieweg, Bonn. 
(1996) 

[3] Kobayashi T., Singular solutions and prolongation of holomorphic solutions to nonlinear 

differential equations, Publ. RIMS, Kyoto Univ., 34(1998), 43-63. 
[4] Kichenassamy S. and Littman W., Blow-up surfaces for nonlinear wave equations, I, Comm. 

Partial Differential Equations, 18(3&4)(1993), 431-452. 
[5] Kichenassamy S. and Littman W., Blow-up surfaces for nonlinear wave equations, II, Comm. 

Partial Differential Equations, 18(11)(1993), 1869-1899. 
[6] Kichenassamy S. and Srinivasan G. K., The structure of WTC expansions and applications, 

J. Phys. A: Math. Gen. , 28(1995), 1977-2004. 
[7] Tahara H., On the singularities of solutions of nonlinear partial differential equations in the 

complex domain, "Microlocal Analysis and Complex Fourier Analysis", World Sci. Publ., 

Hackensack, NJ (2002), 273-283. 
[8] Tahara H., On the singularities of solutions of nonlinear partial differential equations in the 

complex domain, II, "Differential equations & asymptotic theory in mathematical physics", 

Ser. Anal., World Sci. Publ., Hackensack, NJ (2004), 343-354. 
[9] Uchikoshi K., Singular Cauchy problems for quasilinear equations of order two, J. Math. 

Pures Appl. (9) 83(9)(2004), 1151-1178. 
[10] Yamane H., Nonlinear Cauchy problems with small analytic data, Proc. Amer. Math. Soc. 

134(11)(2006), 3353-3361. 

Department of Physics, Kwansei Gakuin University, Gakuen 2-1, Sanda, Hyougo 669- 
1337, Japan 

E-mail address: yamaneOksc.kwansei.ac.jp 



